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Abstract. We report on the design and performance of a highly sensitive combined 3 He/Cs magnetometer 
for the absolute measurement of magnetic fields. The magnetometer relies on the magnetometric detection 
of the free spin precession of nuclear spin polarized 3 He gas by optically pumped cesium magnetometers. 

We plan to deploy this type of combined magnetometer in an experiment searching for a permanent 
electric dipole moment of ultracold neutrons at the Paul Scherrer Institute (Switzerland). A prototype 
magnetometer was built at the University of Fribourg (Switzerland) and tested at Physikalisch-Technische 
Bundesanstalt (Berlin, Germany). We demonstrate that the combined magnetometer allows Cramer-Rao- 
limited field determinations with recording times in the range of 10 ~ 500 s, measurements above 500 s 
being limited by the stability of the applied magnetic field. With a 100 s recording time we were able to 
perform an absolute measurement of a magnetic field of « 1 fiT with a standard uncertainty of AB ~ 60 fT, 
corresponding to AB/B <6xl0~ 8 . 

PACS. 07.55.Ge Magnetometers for magnetic field measurements - 32.30.Dx Magnetic resonance spectra 
- 78.20.Ls Magneto-optical effects - 42.62.Fi Laser spectroscopy 


1 Introduction 


A new experiment searching for a permanent electric 
dipole moment of the neutron (nEDM) is currently being 
developed at Paul Scherrer Institute (PSI), Switzerland 
PQ. In the experiment, the degeneracy of the neutron’s 
magnetic sublevels is lifted by the interaction of the neu¬ 
tron’s magnetic moment = gifijsil/h = Uni/I with 
a static magnetic field B 0 , where I is the neutron’s an¬ 
gular momentum. Additionally, a static electric field £ is 
applied parallel or antiparallel to Bq. In case the neutron 
has an EDM d n = d n I/I , the external field interaction 
Hamiltonian is given by 

H ext = -^IB 0 -j-I- £ - (!) 

The Larmor precession frequencies 
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in transverse (perpendicular to I) magnetic fields are 
given by the sublevel splitting for parallel (ft) and an¬ 
tiparallel (f).) magnetic and electric fields (2|. In the PSI 
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experiment the precession frequency of spin-polarized ul¬ 
tracold neutrons is measured by Ramsey’s method of 
(time-) separated oscillatory fields [5]. A measurement cy¬ 
cle takes « 400 s and consists of filling the neutron stor¬ 
age vessel, running a Ramsey cycle, emptying the vessel 
and measuring the neutrons’ spin polarization. The nEDM 
experiment measures whether or not the (magnetic) spin 
precession frequency is altered by an electric field applied 
along the magnetic held. 

In practice one compares the precession frequencies in 
experiments with (ft) and (f|) configuration. From the 
difference frequency 

Wft — uqq = (£ff + £\T) + -j— (Ho-j-f — ^074) 5 ( 3 ) 

the nEDM is inferred. From Eq. (|3| it is clear that the 
magnetic held has to be precisely controlled and known 
during the precession time of the neutrons (typically 
~ 180 s), hence the need for very accurate magnetome¬ 
ters. An uncorrected statistically huctuating difference 
-B-jq — -B-fi 7 ^ 0 would worsen the statistics for the nEDM 
experiment since the electric-held-induced difference in 
precession frequencies (if any) is many orders of magni¬ 
tude smaller than the Larmor frequency itself. Even more 
severe, a change in the magnitude of the applied mag- 
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netic field that is correlated to the electric field direction 
would be misinterpreted as an nEDM, if not corrected for. 
Moreover, magnetic field gradients may lead to geometri¬ 
cal phase effects j2] that could also mirnick an nEDM. 
A high-sensitivity and high-accuracy measurement of the 
magnetic field and its variation over the neutron stor¬ 
age volume during the free evolution time is therefore of 
crucial importance for controlling and suppressing several 
major systematic errors and ensuring good statistics in 
the nEDM experiment. 

2 Magnetometry in the PSI-EDM experiment 

The current nEDM experiment at psipp deploys two 
types of optically pumped magnetometers for measuring 
the temporal and spatial variations of the magnetic field 
in and around the neutron storage volume, viz., a Hg 
comagnetometer and an array of 16 Cs magnetometers 
(CsOPMs). We first note that the Hg and the Cs mag¬ 
netometers as well as the magnetometer described in this 
work are all scalar magnetometers that measure only the 
modulus of the magnetic field vector at the sensor’s loca¬ 
tion. The working sensitivities of the currently deployed 
systems (Cs, Hg) are both ~ 100 fT in a 100 s measure¬ 
ment time. This is sufficient for the ongoing phase of the 
project since the statistical sensitivity of the nEDM mea¬ 
surement is currently limited by the neutron counting 
statistics. The uncertainties arising from this limitation 
currently allow a tolerance of magnetic field fluctuations 
up to « 100 fT during one Ramsey cycle without loss of 
sensitivity. The Hg co-magnetometer [3] yields a volume- 
averaged value of the field in the neutron bottle, impor¬ 
tant for normalization of the neutron precession frequency. 
The co-magnetometer uses (nuclear) spin-polarized 199 Hg 
vapor that occupies the same storage volume as the neu¬ 
trons, and the magnetic field is inferred from the frequency 
of the free spin precession (FSP) of the Hg’s spin polar¬ 
ization. The Hg-FSP is monitored by recording the (time 
dependent) transmitted power of a resonant circularly po¬ 
larized light beam traversing the Hg vapor. Strictly speak¬ 
ing, the spin precession of the Hg atoms is only ‘quasi’- 
free, since the read-out light beam may affect the spin 
precession frequency (e.g., by the light shift effect m , 
thereby limiting the magnetometer’s accuracy. The rela¬ 
tive systematic shift arising from this effect for the given 
parameters of the current nEDM experiment has been es¬ 
timated to be on the order of 6 x 10 -9 [3]. A limiting 
factor to the accuracy of field measurements with the Hg 
magnetometer is the Hg gyromagnetic ratio, which is only 
known with a relative accuracy of ~ 10 " 6 0 , 0 . 

The Cs magnetometer array builds on an optically de¬ 
tected magnetic resonance process [8], in which the fre¬ 
quency of a weak applied oscillating magnetic field is 
made identical (‘locked’) to the Cs atoms’ Larmor pre¬ 
cession frequency using a feedback loop [5] . The Cs mag¬ 
netometers offer the possibility to access information on 
the spatial field distribution outside the UCN precession 
volume, needed to control magnetic field gradients. The 
Cs magnetometers are also prone to possible light shift 


effects, but their accuracy may be more seriously affected 
by phase errors (and the stability thereof) in their feed¬ 
back electronics. Moreover, because of the hyperfine inter¬ 
action, the Cs magnetometer readings are affected by the 
quadratic Zeeman effect and, being driven magnetome¬ 
ters (in contrast to magnetometers based on free preces¬ 
sion), their interaction with the magnetic resonance driv¬ 
ing rf field introduces a systematic frequency shift (Bloch- 
Siegert shift [TO]). The systematic effects mentioned above 
lead to sensor-specific offsets of the magnetometric read¬ 
ings that are generally unknown, not necessarily constant 
and may depend on other experimental parameters. This 
spoils the direct comparison of the absolute field value 
given by different sensors and limits the ability to obtain 
a consistent picture of the magnetic field inside the appa¬ 
ratus necessary to suppress systematic errors in the EDM 
measurement such as those mentioned in Sec. [l] caused by 
the geometric phase effects [2j. 

The combined 3 He/Cs magnetometer described hereafter 
offers an important complement to the magnetometers dis¬ 
cussed above. It is based on recording the FSP of nuclear 
spin polarized 3 He gas, by detecting the time dependent 
magnetic field produced by the processing (and decay¬ 
ing) 3 He magnetization with an arrangement of several 
Cs magnetometers. 3 He FSP detection by external mag¬ 
netometers provides an indirect optical readout nonper- 
turbative to the 3 He FSP, avoiding possible systematic 
effects as they may occur by the read-out beams or the 
feedback control in the Hg and Cs magnetometers, re¬ 
spectively. In the absence of magnetic field gradients the 
Larmor precession frequency of 3 He is thus an absolute 
measure of the magnetic field inside the magnetometer 
cell. When gradients are present, the measured precession 
frequency corresponds to the volume averaged field in the 
cell. The details of this averaging process will depend on 
the dynamic regime in the cell m- On the other hand, 
the systematic effects affecting the CsOPMs are irrelevant 
for their use as readouts for 3 He FSP. The fact that the 
3 He gyromagnetic ratio is known with a relative precision 
of 2.5 xlO -8 [7] makes 3 He a promising candidate as a ref¬ 
erence for magnetic field measurements, as suggested in 
(T21 . One application that is envisioned for the next stage 
of the nEDM experiment (n2EDM) is a 3 He- “quasi comag¬ 
netometer”. Since in this experiment it is not possible to 
have 3 He cohabiting inside the neutron bottle, flat cylin¬ 
drical magnetometer vessels will be installed above and 
below the cylindrical neutron precession chamber. These 
vessels will have the same geometrical cross section as 
the precession chamber and will thus be traversed to first 
order by the same magnetic flux. An array of CsOPMs 
around the vessels will be used to detect the 3 He FSP 
signal. The magnetic field measurement from the 3 He/Cs 
magnetometer, performed simultaneously with the nEDM 
measurement proper, can then be used, along with the Hg 
co-magnetometer, to normalize the neutron precession fre¬ 
quency and to correct for field changes |13j.|14|. From the 
double-chamber type of geometry, additional information 
on magnetic field gradients can be obtained. Another pos¬ 
sible application is an array of several compact 3 He/Cs 
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magnetometers to measure the magnetic field at differ¬ 
ent spatial positions around the precession chamber. Since 
3 He/Cs magnetometers do not suffer from the systemat- 
ics discussed above, and offsets which they may induce, 
readings from different 3 He/Cs sensors can be more easily 
compared and their use can thus improve field stabiliza¬ 
tion and gradient control in the experiment. In this paper, 
we describe a prototype 3 He/Cs magnetometer based on a 
small, sealed spherical 3 He cell surrounded by Cs magne¬ 
tometers. The prototype is not adapted to later use in the 
nEDM experiment, but designed as a versatile test device 
for investigation of the combined magnetometer concept. 

2.1 Preparation and detection of polarized 3 He 

3 He gas in a spherical glass cell is polarized by optical 
pumping of the helium atoms in the metastable 2 3 Si state 
populated by collisions in a weak glow discharge, using 
circularly polarized 1083 nm light from a 2W ytterbium- 
doped fiber laser. The electronic spin polarization is trans¬ 
ferred to nuclear spin polarization of the groundstate 
atoms by metastable exchange. Metastable exchange opti¬ 
cal pumping (MEOP) is a method that is well studied and 
described in detail, e.g., in [15] and |16j. Conversely to Cs 
and Hg magnetometers, the spin polarization (and hence 
its dynamics) cannot be detected by optical means be¬ 
cause of the large optical excitation energy (20 eV) of the 
3 He ground state for which no (convenient) light sources 
are available. 

The rotating magnetization mne that is associated 
with the precession of the nuclear spin polarization can 
be detected by pick-up coils which are very ineffi¬ 
cient at the low precession frequency (~ 30 Hz) occuring 
here. 

SQUIDs (superconducting quantum interference de¬ 
vices) are highly sensitive magnetometers that have been 
deployed for detecting the 3 He FSP m- However, the 
additional technical complexity associated with the cryo¬ 
genic cooling needed for the SQUID operation is an ob¬ 
stacle for operating these magnetometers under the ex¬ 
perimental conditions of a room temperature nEDM ex¬ 
periment. In 1967, Cohen-Tannoudji and co-workers m 
have demonstrated the suitability of a discharge lamp 
pumped alkali (Rb) vapor magnetometer for detecting the 
FSP of 3 He nuclei. Based on the results reported in that 
paper, we estimate their magnetometric sensitivity to be 
« 80pT • s 3 / 2 /Tm 2 , assuming a Cramer-Rao limited per¬ 
formance (c.f. Sec. [5]) in a measurement time T m . 

2.2 3 He/Cs magnetometer principle 

During optical pumping the 3 He cell is exposed to a ho¬ 
mogeneous magnetic field B 0 oriented along the pump 
laser beam. The oriented nuclear magnetic moments give 
rise to a macroscopic magnetization mn e that produces 
a magnetic dipole-like field B n e outside of the cell. One 
readily estimates that the field from a 100 % polarized gas 
at 1 mbar is on the order of 200 pT on the outside surface 


of the cell. After optical pumping, the laser and the gas 
discharge are turned off and a 7t/2 rf-pulse is applied to 
the cell in order to flip mn e to a plane perpendicular to 
B 0 , upon which mjj e starts freely precessing around Bq 
at the 3 He Larmor frequency 

^He = 7He |-Bq| , (4) 

where 7 He /27r = 32.43410084(81) Hz///T 
(Ayne/7He =2.5xl0 -8 ) is the gyromagnetic ratio of the 
3 He nucleus J7|. The precessing (and decaying) magnetiza¬ 
tion produces at the position r c s of each Cs magnetometer 
a magnetic field i?He(^CsU) with time-dependent ampli¬ 
tude and orientation. 

The Cs magnetometers are scalar magnetome¬ 
ters, i.e., they measure the modulus B(r,t) = 
|-B 0 (r,f) + -BHe(u t)\ of the total field at their location. 
Since and Bq is nominally constant in time one 

has B ss Bq + B 0 • £?He(t), so that the CsOPMs are, 
to first order, only sensitive to the component SB X of 
the 3 He-FSP field along the applied magnetic field Bq. A 
simple calculation shows that—for a given distance rc s — 
this time dependent projection has a maximum amplitude 
when the sensors are located on a double cone with a 
half-opening angle of y>=45° with respect to Bq. In the 
prototype described below, the centers of all 8 CsOPM 
cells were located on that double cone, and the relative 
azimuthal positions on the cones determine the phase re¬ 
lations between the individual FSP signals detected by the 
different CsOPMs (compare Fig. IT]). By pairwise subtrac¬ 
tion of CsOPM signals that are dephased by n, common 
mode magnetic noise components (such as magnetic fields 
oscillating at the 50 Hz line frequency) that are in-phase 
on both sensors can be strongly suppressed in the differen¬ 
tial signal, while increasing the signal of interest, as shown 
in Sec. |5.2| For the chosen geometry of the combined mag¬ 
netometer, 8 such gradiometer pairs can be formed. 

2.3 The Cs magnetometers 

The CsOPMs used in this study are laser pumped double¬ 
resonance magnetometers, operated in the M x configura¬ 
tion [8. A sketch of a CsOPM is shown in Fig. [2] (left). 
Light from a Toptica DLPro diode laser with a frequency 
actively stabilized to the F g =4—>F e =3 transition of the 
Cs-Di line at 894 nm is delivered to each magnetometer 
module via a multimode fiber. The light from each fiber is 
collimated and given a circular polarization, after which 
it passes through a room-temperature paraffin-coated [20] 
30 mm diameter Cs vapour cell. The cell coating ensures 
a long-lived coherence (~ 30 ms) of the spin polarization 
created by optical pumping. The light exiting the cell is 
detected by a photodiode (PD) that measures the Cs va¬ 
por’s optical transmission. The propagation direction kc s 
of the incident light makes an angle of 45° with respect to 
Bq for all CsOPMs since this yields maximal sensitivity 
[8]. A weak magnetic field (rf field) parallel to kc s oscil¬ 
lating at frequency w r f resonantly drives the precession of 
the Cs vapour’s magnetization, thereby modulating the 
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Fig. 1. 3 He spin sample cell with its magnetic dipole field in 
the (a;,y)-plane. The direction of the applied Bo defines the 
x axis, so that the 3 He magnetization precesses in the ( y,z ) 
plane. The Cs sensors shown at their different positions record 
the FSP signal phase-shifted by n. The 45° cones of highest 
sensitivity are denoted by the dashed lines. 


vapour’s absorption coefficient and hence the photodiode 
signal [Sj. The (transimpedance-) amplified PD signal is 
demodulated by a dual channel digital lock-in amplifier 
referenced to w r f. Figure [3] shows the dependence of the 
amplitude and phase of a CsOPM sensor on the detuning 
from the resonance S = w r f — uqs ■ The resonance occurs 
at the Cs Larmor frequency w r f = ujqs = 7 Cs | -Bo | , where 
7Cs ~ 27 t x 3.5kHz//rT is the Cs ground state’s gyromag- 
netic ratio. In the currently employed mode of operation 
in the nEDM experiment the rf frequency is tuned near 
the line center where the phase has a linear dependence 
on the frequency detuning and the phase signal is used to 
drive a voltage-controlled oscillator generating the oscilla¬ 
tory voltage for the rf coils. This mode of operation repre¬ 
sents a feedback loop that keeps w r f locked to the Cs Lar¬ 
mor frequency 0. When multiple CsOPMs are operated 
in close spatial vicinity one sensor may be parasitically 
driven by the rf of a neighboring magnetometer. To avoid 
this effect known as cross talk, it is advisable to drive all 
CsOPMS at the same common frequency. In the measure¬ 
ments described here the CsOPMs were thus driven at a 
single, constant frequency u > T f close to uic s and the photo¬ 
diode signals were demodulated at that fixed rf-frequency 
using digital lock-in amplifiers. The oscillatory magnetic 
field Bn e (t) then leads to an oscillation of the CsOPM’s 
phase signal at wne, with an amplitude proportional to 
i?He, as long as L?He7 cs^A*;, where Auj is the linewidth 
of the resonance of Fig. [3j We note that the proportional¬ 
ity factor between the phase response of the CsOPM and 
i?He depends on the Cs cell properties and the parameters 


V 


rf coils photodiode 



lens A/4 

polarizer Cs cell A 6, 



Fig. 2. (Color online) Schematic drawing of the CsOPM (left) 
and CAD drawing of the combined 3 He/Cs magnetometer 
(right). The spherical 3 He cell (magenta) in the middle of 
the cubic structure is surrounded by eight CsOPMs (blue) on 
the edges of the cube, in which the rf coils are laid out on 
printed circuit boards (yellow). The total dimensions of the 
combined magnetometer are ~ (10cm) 3 . One corner cube and 
two CsOPMs are left out for better visibility. 


of operation. For the CsOPMs driven at fixed frequency, 
bandwidth limitations arising from the lifetime of the Cs 
polarization and the lock-in demodulator filter have to be 
considered. The scaling of the phase signal to magnetic 
units which we will do in the following to make the re¬ 
sults more intuitively accessible does not correct for these 
effects, e.g., the reported amplitudes do not reflect the 
true values of Bn e . The true FSP amplitudes are actually 
roughly a factor of five larger. 


3 The prototype magnetometer 

A prototype of a combined 3 He/Cs magnetometer for 
studying the 3 He FSP readout with laser pumped 
CsOPMs was built in Fribourg. It consists of a 70 mm di¬ 
ameter spherical cell filled with 3 He at a pressure of 1 rnbar 
that is fixed in the center of a mechanical structure hold¬ 
ing eight laser-pumped CsOPMs mounted symmetrically 
on the edges of a cube as shown in Fig. [2] (right), thereby 
fulfilling the optimal sensitivity criterion (45° cone) dis¬ 
cussed above. The distance between the 3 He cell center 
and the Cs cell centers is 50.5(5) mm. The He cell car¬ 
ries two electrodes, each consisting of a spiral of copper 
foil glued to the outside of the cell that are driven by a 
1.2 MHz sinusoidal voltage (amplified by a Tesla trans¬ 
former) to ignite and sustain a weak gas discharge in the 
cell. The electrodes were designed to achieve a homoge¬ 
neous illumination of the cell volume by the gas discharge 
while allowing optical access to a large part of the cell’s 
surface to permit the 3 He pump laser beam to traverse 
the cell. 

Figure[2](left) shows details of a single CsOPM sensor. 
Each sensor carries its own pair of rf coils and the coils 
of all 8 sensors were driven at the same, constant rf fre¬ 
quency during measurements. This mode of operation is 
only possible when the magnetic held gradients over the 
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Fig. 3. Amplitude (top) and phase (bottom) response of a 
CsOPM when sweeping the rf frequency lj t {. The dashed lines 
represent the change of the CsOPMs Larmor frequency u>cs 
and phase response corresponding to a 1 nT variation of the 
magnetic field Bq. 
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whole structure are sufficiently small so that the individual 
Larmor frequencies of all 8 sensors differ by amounts that 
are much less than the Cs magnetic resonance linewidth 
Alu. During the measurements described here, the local 
Larmor frequencies of two CsOPMs i,j £ {1,8} differed 
by less than (w* — Wj)/(2n) < 1.8 Hz. 

The apparatus depicted in Fig. [2] (right) is surrounded 
by large coils (Helmholtz configuration, ~ 30 cm diame¬ 
ter, not shown in the figure) for producing magnetic fields 
perpendicular to B 0 . They were used to start the FSP 
by flipping the 3 He magnetization by resonant rf-pulses 
following the pumping process. Having in mind the later 
mounting of the device in the vacuum chamber of the 
nEDM spectrometer at PSI, all components were man¬ 
ufactured from nonmagnetic and vacuum compatible (low 
outgassing) materials. 


4 Measurements 


A key issue of the present study was the determination 
of the intrinsic magnetometric sensitivity of the combined 
3 He/Cs magnetometer. Since the stability of the applied 
magnetic held sets a limit on the ability to determine the 
magnetometer sensitivity (c.f. Sec. 5.4 for details), that 
held has to be kept as stable as possible. After initial 
tests in Fribourg the measurements reported below were 


carried out in the magnetically shielded room BMSR-2 
m at the Physikalisch-Technische Bundesanstalt (PTB) 
in Berlin, Germany. The BMSR-2 room is one of the mag¬ 
netically most quiet and stable places on earth. It consists 
of a 7 layer MU-metal magnetic shield and an additional 
aluminum layer. It features a built-in multi-channel sys¬ 
tem of (vector) SQUID magnetometers that were operated 
together with our test equipment. A three axis Helmholtz 
coil system [22] was used to produce the B 0 held inside 
BMSR-2. The coil was driven by a commercial low-noise 
current source (Magnicon, CSE-1) delivering a current 
of 19 mA yielding a homogeneous magnetic held |i?o| of 
ss 1 fj,T in the center of the coil. The 3 He optical pumping 
light from the 1083 nm laser was brought into the chamber 
by an optical hber, after which the beam was expanded 
by a telescope in order to illuminate the whole accessible 
cross section of the 3 He cell. The light was circularly polar¬ 
ized with a polarizing beamsplitter cube followed by a A/4 
plate. After traversing the 3 He cell the pump beam was 
back-reflected by a mirror for a second passage through 
the cell, thus increasing the pumping efficiency. The fluo¬ 
rescence from a 3 He reference cell located outside BMSR-2 
was used to monitor and manually adjust the laser wave¬ 
length during optical pumping. 

The 894 nm light for operating the 8 CsOPMs was de¬ 
livered to the sensors by eight multimode fibers. A board 
containing 8 transimpedance amplifiers (mounted inside 
BMSR-2) pre-amplffied the photodiode currents, and the 
ensuing voltage signals were transmitted to the data ac¬ 
quisition (DAQ) system located outside of the chamber. 
The raw signals as well as the demodulated PD signals 
from a set of six CsOPMs were recorded simultaneously by 
(nominally) identical DAQ channels. The remaining two 
CsOPM signals were recorded by a separate DAQ system 
and were not used for the analysis presented here. The 
timebase of the DAQ system was referenced to a rubid¬ 
ium atomic clock (SRS PRS10). 

Figure [4] shows the demodulated phase signals of six 
CsOPMs after rescaling to magnetic field units. The ~ 
36 Hz oscillation from the 3 He FSP is clearly visible with 
no additional filtering applied to the phase data. If not 
specified otherwise, the 36 Hz precession signals from the 
3 He atoms will be referred to as FSP signals in what fol¬ 
lows. The FSP signals in Fig. [4] are clearly of varying 
quality. While Csl (upper left in Fig. El) exhibits a rel¬ 
atively good signal to noise ratio (SNRj, Cs8 (lower left 
in Fig. 0 is obviously performing much worse. These dif¬ 
ferences are due to different Cs-cell qualities and CsOPM 
drive parameters. A criterion to quantify the sensors in¬ 
trinsic sensitivity limitation taking cell quality and drive 
parameters into account is the noise equivalent magnetic 
field (NEM) (c.f. Sec. [6] for details). While for Csl we find 
a quite good value of NEMc s i « 34fT/\/Hz, Cs8 yields 
NEMc S 8 « 158 fT/-s/Hz. For the measurements presented 
here, Csl was always performing significantly better than 
the other CsOPMs. While the drive parameters can be op¬ 
timized during a measurement, preselection of high qual¬ 
ity Cs cells for the combined magnetometer is of crucial 
importance. 
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Fig. 4. Raw phase signals of the six CsOPMs used in the 
analysis, rescaled to magnetic units. All subplots have identical 
amplitude and time scales. The ~ 36 Hz oscillation originating 
from the 3 He FSP is clearly visible. The amplitudes and signal 
to noise ratios differ due to individual bandwidth limitations 
of the individual CsOPMs. Note that the amplitudes do not 
reflect the true magnitude of Bu e - 


5 Data analysis 

The data were analyzed off-line using dedicated Mathe- 
matica [23; codes. 

5.1 Relaxation time 

As discussed in EL the transverse spin relaxation time 
T 2 of nuclear spin polarized 3 He atoms is strongly affected 
by the presence of magnetic field gradients. The very long 
T 2 time that can be achieved in high quality glass cells is, 
in general, limited by field inhomogeneities. We recorded 
3 He FSP signals for time periods of slightly more than 10 
hours. The data, recorded at a sampling rate of 450 Hz, 
was split into 44 s long segments over which the FSP am¬ 
plitude can be assumed to be constant. Each subset was 
then individually fitted using a sinusoidal function 

fit) = a 0 ff + aft) sin (uj He t + <t >), (5) 

where a 0 ff is the offset field at the individual sensor’s po¬ 
sition and a(t) the constant FSP amplitude for the re¬ 
spective subset. The time dependence of these amplitudes 
for a single CsOPM is shown in a semi-logarithmic plot 
in Fig. [5] that illustrates the exponential character of the 
decay. The FSP signal is described by 

s(t) = ao e~ t i T2 sin(wf/ e t + <j>). (6) 

From a fit of the data shown in Fig. [ 5 ] (from a single 
CsOPM) by the function 

a(t) = a 0 e~ t/T2 (7) 

we infer a decay time T 2 of 13173(4) s ~ 3.6 h. Closer in¬ 
spection of the fit residuals reveals a small imperfection of 
the fit (red curves, middle graph of Fig. [5]) . This can be ex¬ 
plained by variations of magnetic held gradients induced 
by the magnetic relaxation of the /Lt-metal shield (essen¬ 
tially the innermost layer) after closure at the beginning 
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Fig. 5. (Color online) Measurement of 3 Hc T 2 time. Fit of 
single exponentially decaying function Eq. ([7|l and combined 
function Eq. j8j, taking into account the relaxation of magnetic 
field gradients (top). The fit residuals are shown for the simple 
function (middle, red) and the more complex one (bottom, 
green). Note that the amplitudes -as in Fig. [4] do not reflect 
the true magnitude of Bn e - 


of the measurement. This effect has already been observed 
in BMSR-2 and was described in [23]. An alternative fit 
function, 


a(t) = ai e t/T i 1) +a 2 e tlT ? ) (8) 

empirically takes into account the shield relaxation by in¬ 
troducing a second time constant . The fit of Eq. ([8]) 
yields Gaussian-distributed residuals and a decay time 
X ^ of 13505(23) s (green curve, lower in Fig. [ 5 ]). The same 
analysis, performed with all six simultaneously running 
CsOPMs yields the decay rates visualized in in Fig. [6] 
As expected, the values for (and respectively) 
measured by different CsOPMs agree within their un¬ 
certainties. The mean decay times of the 3 He polariza¬ 
tion, calculated from all six CsOPM measurements, are 
= 13532(17) s and = 6621(183) s respectively. 
Although much longer decay times have been reported in 
the literature (see, e.g., 125 ]), the value achieved here is 
largely sufficient in the context of the present study. 

In a constant magnetic field the 3 He FSP can be repre¬ 
sented by a decaying single tone oscillation as introduced 
in Eq. ©■ Assuming that the data have only white Gaus¬ 
sian noise Q and neglecting the shield relaxation described 
above, we can model the experimental signal as a discrete 
time series of equi-spaced data points 

S n = ae nT i T2 sin(wnT + 4> 0 ) + Q(n ), (9) 
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CsOPM 



Fig. 6. (Color online) Decay times of 3 He spin polarization 
simultaneously measured by multiple CsOPMs. The red hor¬ 
izontal lines denote the weighted mean and the shaded band 
around it gives the la confidence region. Compare Eq. |8| for 
meaning of variables. 



0.1 0.2 0.5 1 2 5 10 20 50 100 


r 


Fig. 7. Plot of CRLB degradation function C(r) from Eq. (131 
for a damped sine wave as a function of r = Tij Tm , the ratio 
of decay to measurement time. 


where T = (f sr ) _1 is the inverse of the sampling rate f sr , 
i.e., the spacing between consecutive points in the time se¬ 
ries and S(n) the Gaussian noise contribution to the n-th 
data point. The noise is completely characterized by its 
power spectral density p 2 or variance <jg = p 2 fbw where 
fbw = fsr/ 2 is the bandwidth of the measurement. The 
precision of the frequency determination of such a coher¬ 
ent signal over a given measurement time Tm is funda¬ 
mentally limited as described by information theory. It 
has first been studied by Cramer and Rao [263, El who 
derived a lower bound for the frequency estimation vari¬ 
ance of a signal of constant amplitude (CRLB) [55] ■ The 
corresponding bound for a damped oscillation was derived 
in [15] and reads 


of > 


(2t r) 2 SNR 2 T 3 , 


■C(T m ,T,T 2 ), 


( 10 ) 


with 


C(T m , T, T 2 ) = 

Tl _ (1 ~ e ~ 2r / T2 ) 3 (1 — a ) 

6 T 3 e _2T / T ' 2 (1 - a) 2 - (T m /T) 2 a(l - e - 2T / T2 ) 2 

( 11 ) 


where 


a = e ~ 2TM/T2 


( 12 ) 


InEq. (101, SNR = a r 


ls /p represents the signal to noise- 
density ratio with a rms = a/yj 2 being the rms-amplitude 
of the FSP signal. C(Tm,T, T 2 ) is a factor that takes 
damping into account. For sufficiently high sa mplin g rates 
(T <C 27 t/o iHe), an approximative form of Eq. (11) that is 
independent of T can be found and reads 


C(r) = 


s 2 / r - 1 


3r 3 cosh ( 2 ) — 3r (r 2 + 2) ’ 


(13) 


where r = T 2 /Tm is the ratio of the decay and measure¬ 
ment time. A plot of this function is shown in Fig. [7] For 


vanishing damping, T 2 Tm the factor C(Tm,T,T 2 ) be¬ 
comes unity, thus reproducing the result of [2S|. Since the 
decay time T 2 in our system is ~ 13,000 s, one sees from 
Fig.0 that C is «1 for measurement times up to several 
thousand seconds and can thus be neglected. It turns out 
that up to integration times of several hundreds of seconds 
the fit results do not differ when the decay is neglected. 
This justifies the use of the simple fit function ([ 5 ]) up to 
relatively long integration times for which the signal am¬ 
plitude can be considered constant. 


5.2 Common noise suppression 

Noise suppression effects in differential CsOPM signals as 
described in Sec. |2.1| were investigated. For this, the phase 
signal from one CsOPMs (Cs5) was subtracted from a sec¬ 
ond CsOPM signal that is, by construction, dephased by 7 r 
(Cs4). Both signals carry perturbations from the 50 Hz line 
frequency, as shown in the Fourier spectra in the top row 
of Fig. [ 8 ] These perturbations are in-phase on both signals, 
whereas the 3 He FSP signal is dephased by 7 r, as visible 
in the lower left plot of Fig. [ 8 ] In the differential signal, 
the 50 Hz perturbation has vanished, as evidenced by the 
Fourier spectrum on the lower right of Fig. [ 8 ] It is expected 
that the (random) noise amplitude spectral densities of 
the two signals p 4 = 48fT/\/0z and p$ = 59fT/-\/Hz add 
quadratically in the differential signal pdiff = \Jp\ + p\ > 
assuming no correlation between the white noise contri¬ 
bution to both signals. One also expects the signal ampli¬ 
tudes 04 = 3.58pT rms and at 5 = 3.86pT rms to add in the 
combined signal, Odiff = 04 + 05 . The expected signal to 
noise ratio of the differential signal can thus be written as 

SNRd iff = , (14) 

VpI + Pi 

yielding SNRdiff ~ 98\/Hz for the values given above. 
Analysis of a Fourier spectrum of the differential signal 
shows that only a reduced SNR of ~ 85 \/Hz is observed 
in the signal. A closer analysis reveals that the amplitudes 
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frequency (Hz) 



frequency (Hz) 



frequency (Hz) 


Fig. 8. Upper row: Fourier spectra of two CsOPMs dephased by ir. Perturbations from the 50 Hz line frequency and around 
~ 22Hz are visible. Lower row: Time series of both signals, dephasing is visible (lower left), Fourier spectrum of differential 
signal. The 50 Hz pertubation has vanished but the ~ 22Hz perturbations persist and are even increased (lower right). Each 
Fourier spectrum contains 700 s of measurement data. 


add, as expected, to ~ 7.43 pT, while the noise of the dif¬ 
ferential signal is increased more than expected to a level 
of « 87 fT/v / Hz. This discrepancy can be explained by 
non-Gaussian perturbations which are correlated between 
the two signals. The presence of such noise components 
is witnessed by the double-peak around ~ 22 Hz in the 
Fourier spectra of Fig. [8] Such perturbations may origi¬ 
nate from higher order magnetic field fluctuations (gradi¬ 
ent oscillations) that affect the individual sensors located 
in different spatial positions differently, but in a correlated 
manner. Such noise processes are not suppressed by the 
gradiometer, their amplitudes might even add up in the 
combined signal, as visible in Fig. [8] (lower right). Never¬ 
theless, the suppression of the common noise component 
was successfully demonstrated. This technique becomes 
an important tool when in-phase perturbations are strong 
and may lead to systematic errors in the extraction of the 
3 He Larnror frequency. 


5.3 Magnetic field measurements 


the BMSR-2 chamber by analyzing consecutive 100 sec¬ 
ond long time series of the continuously recorded phase 
signal. We analyzed data from the six CsOPMs that have 
recorded simultaneously the same 3 He FSP. The average 
frequency of each time series and its standard error were 
extracted by fitting Eq. ^5]) to the data. Special care was 
taken to ensure that the nt routine correctly estimates the 
standard errors. Since the phase data undergoes filtering 
in the lock-in amplifier, the noise might not be purely 
Gaussian anymore. To prevent the fit routine from under¬ 
estimating the error, the variance was extracted from a 
Fourier spectrum of each dataset and explicitly imposed 
on the fit. The magnetic field Bi and its uncertainty ABi 
were calculated for each CsOPM using Q. The weighted 
average B and its uncertainty AB were calculated for the 
readings of all individual CsOPMs according to 


B = 


E 


AB 2 


E 


AB 1 


with standard error 


(15) 


In order to demonstrate the performance of the com¬ 
bined magnetometer we measured the magnetic held in 




AB ? 


- 1/2 


(16) 
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Fig. 9. (Color online) Time evolution of the magnetic field 
in the BMSR-2 chamber at PTB, estimated from consecutive 
100 s subsets from data segments of six CsOPMs simultane¬ 
ously detecting the same 3 He FSP. Each plot point displays 
an individual CsOPM reading Bi and its uncertainity ABi. 
The common ordinate is shifted by ~ 1.1 /rT. The solid line 
(red) and the shaded band represent B and its corresponding 
la confidence values B ± AB. 


Equations (15) and (16) are only valid assuming a con¬ 
stant magnetic field which is justified for 100 s subsets 
of data as shown in Sec. |5.4[ The uncertainities ABi of 
the field estimations from data of an individual CsOPM 
range between 100 and 350 fT in 100 seconds, depending 
on the signal to noise ratio of the respective CsOPM and 
the field fluctuations during the measurements. The error 
AB of the weighted mean is typically well below 70 fT for 
these 100 second time slices. 

Figure [9] shows, that the magnetic field in the BMSR-2 
chamber drifts by several pT on a time scale of « 1000 s. 
We will see below that these instabilities limited the de¬ 
termination of the combined magnetometer’s sensitivity. 


5.4 Experimental determination of sensitivity 

The noise of any magnetometer signal will contain in gen¬ 
eral contributions from both technical and intrinsic (fun¬ 
damental) noise sources. In the best possible case the tech¬ 
nical noise is due to the instability of the applied mag¬ 
netic field. The two sources of noise cannot be distin¬ 
guished in general. In this study we have operated simul¬ 
taneously four different types of magnetometers, viz., the 
3 He read out by CsOPMs, the 3 He read out by SQUIDs, 
the CsOPMs alone, and the SQUID(s) alone. We can make 
use of this fact to distinguish noise processes that are com¬ 
mon to all sensors from processes that affect the individ¬ 
ual sensors. When different types of magnetometers are 
operated in the same magnetic field, the technical noise 
contributions resulting from field fluctuations should be 
correlated, except for instabilities of the magnetic field 
gradients, since the different magnetometers are located at 
different spatial positions. The Allan standard deviation 
(ASD) i29 is a powerful tool to examine noise processes 
and signal stability. A detailed account of its application 



Fig. 10. (Color online) Allan standard deviation of simultane¬ 
ous field measurements with the 3 He/Cs (green), the CsOPM 
(red), the SQUID (blue) magnetometer, and the 3 He/SQUID 
(orange) respectively. For long integration times the sensitiv¬ 
ity of all magnetometers is limited by a common noise process. 
The oscillatory component of the SQUID data is most proba¬ 
bly due to mechanical vibrations of the setup. Note that the 
3 Hc/Cs curve lies deeper than the 3 He/SQUID only due to 
the lower SNR of 3 He/SQUID which is caused by the larger 
distance to the 3 He cell. 


to magnetometric measurements is given by Groeger et 

al- ffl- 

We have extracted the magnetic field from the 3 He / 
SQUID data using the same procedure as for the 3 He/Cs 
signal analysis, i.e., by fitting Eq. ([ 5 ]) to the data. The 
magnetometric field readings from theSQUIDs proper and 
from the Cs magnetometers proper were retrieved from 
the original data using a digital low-pass filter (bandwidth 
~ 25 Hz) that removes the modulation from the 3 He FSP. 
We note that these filtered magnetometer signals do not 
represent absolute field readings since they may be af¬ 
fected by detector specific offsets. Nevertheless, as long as 
those potential offsets are constant, they will not affect 
the ASD of the signals and can thus be used for com¬ 
paring the signal fluctuations. The ASD of the measured 
magnetic field was calculated over a wide range of inte¬ 
gration times for all four types of magnetometers. Typical 
results are shown in Fig. [To] The region of CRLB lim¬ 
ited measurements up to integration times Tm ~ 500 s is 
restricted by the drift of the magnetic field witnessed in 
Fig-i One sees that the ASDs of all sensors, independent 
of the sensor type, end up on the same rising slope for suf¬ 
ficiently long integration times Tm > 500 s. The fact that 
the different ASDs do not overlap perfectly for long inte¬ 
gration times can be explained by fluctuations of magnetic 
field gradients, since the different sensors were not located 
at exactly the same spot. These observations support the 
assumption that the long-term stability of the magnetic 
field limits indeed the determination of the magnetomet¬ 
ric sensitivity. Note that for the data which entered the 
analysis displayed in Fig. 10 the initial 3 He polarization 
was ~ 20 % smaller than for the measurements presented 
in Sec. 15.31 
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In order to experimentally determine the sensitivity 
of the 3 He/Cs magnetometer we selected a measurement 
with a large 3 He polarization and calculated the ASD 
of the magnetic field determined from n>H e following the 
same procedure as above. The result is shown in Fig. m 
the shaded band giving the lcr confidence region for each 
measurement. Where the field estimation process is CRLB 
limited, we find the characteristic r“ 3 / 2 slope of the ASD 
curve, e.g., up to integration times of ~ 55 s. We deter¬ 
mine the experimental sensitivity in this range by fitting 
a function 

_meas meas IT^I 2 / -i n\ 

a B Vb / J M > DC 

to the data (up to r < 55 s). The fitted function is denoted 
by the dashed line in Fig. El The uncertainties of the 
ASD points are used as weights in this fit. We find for the 
sensitivity parameter r)g eas = 107.0(5) pT • s 3 / 2 . 


of the CsOPMs’ photodiode currents and the shotnoise 
of the .B 0 -coilcurrent. We refer to this fundamental noise 
level as the shotnoise limit in the following. The funda¬ 
mental intrinsic noise floor of the measurement will thus 
have two contributions, one from the magnetic field noise, 
Afield, and one from the CsOPM noise, acs ■ Since the two 
Gaussian sources of noise are uncorrelated one can set 

a Q = \! a Gs + Afield ■ ( 18 ) 

As described in [20], the intrinsic noise of a CsOPM de¬ 
pends on the laser power, rf-power and the cell properties 
and can be expressed in terms of the noise-equivalent mag¬ 
netic field (NEM). For optimized rf- and laser power one 
obtains a minimal NEM min of « 12fT/\/Hz for the type 
of cells used. To estimate the NEM in the shotnoise limit 
under given experimental conditions, e.g., non-optimal rf 
and light power, we start from the square-root power spec¬ 
tral density of the photocurrent shotnoise 

Pi,psn = \Z^e.I D c , (19) 

where /dc is the measured DC-photocurrent and e the 
elementary charge. The voltage shotnoise is then given by 

Pv.psn = Pi,psn 9 = \JleIoc 9 1 (20) 

where g = 2.53 x 10 7 V/A is the gain of the used tran¬ 
simpedance amplifier. The phase change 5(f> corresponding 
to a voltage change SU of the input signal is m 



n SU 

dcp = -, 

acs 


where ac s denotes the amplitude of the signal after the 
transimpedance amplifier, and the NEM can finally be 
calculated using 


CO 

< 


10 20 40 60 80 100 200 300 

T (S) 

Fig. 11. (Color online) ASD of magnetic field calculated from 
uiHe measured by the 3 He/Cs magnetometer. The gray region 
around the curve denotes the lcr confidence band for each ASD 
point. A fit of Eq. ( |17| ) to the data up to r < 55 s is displayed 
as the dashed line. 


We note that due to the higher 3 He polarization com¬ 
pared to the data shown in Fig. [lO] the ASD reaches its 
minimum imposed by the instability of the magnetic field 
already at shorter integration times. 

6 Sensitivity in shotnoise limit 

In this section we will quantify the intrinsic noise level that 
is inherent to the measurement process with the given ex¬ 
perimental parameters. We can then compare this value 
to the observed noise floor to judge the quality of our mea¬ 
surement, allowing us to quantify the impact of technical 
noise (such as electromagnetic pick-up, amplifier noise,...) 
on the signals. We will furthermore estimate the ultimate 
achievable sensitivity of the combined magnetometer un¬ 
der optimized experimental parameters. 

The discrete nature of the charges constituting the 
relevant electric currents in the experiment is a funda¬ 
mental source of noise. This shotnoise is assumed to be 
a zero-average Gaussian process. We calculate the intrin¬ 
sic noise level of the measurement due to the shotnoise 


NEM = pcs = — , (22) 

7 Cs 

with r 2 being the transverse relaxation rate of the Cs 
spin polarization. For the best-performing CsOPM under 
true experimental conditions (Csl: I dc = 3.25 pA, ac s = 
1.3 V, r 2 /27T = 6.1Hz) this calculation yields 

NEMc sl ~ 34fT/\/Hz. The photocurrents lr>c were 
repeatedly measured for each sensor during the experi¬ 
ment. The values of ac s and r 2 /2 can be obtained by fits 
to the frequency sweep-responses of the individual sensors 
shown in Fig. [3] which were also repeatedly recorded. 

The second source of noise in Eq. flSj ) comes from 
Gaussian fluctuations of the magnetic fielcD B 0 due to the 
concurrent’s shotnoise. In the shotnoise limit pfi e id is ob¬ 
tained for a given coilcurrent (Ic = 19 mA) and coil con¬ 
stant (gcoii = 60 pT /A) is given by 

p field = \/2 el c g coi i = 4.7 fT/\/Hz . (23) 

Inserting these values into Eq. © leads to the shot¬ 
noise limit under the given experimental conditions of 
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Pg N = 35fT/\/Hz. It is obvious that pg N is dominated 
by Pcs; the contribution from the coilcurrent pg e id play¬ 
ing only a negligible role. This value can be compared to 
the experimentally observed noise level p e g P = 27 fT/\/Hz. 
In order to make this comparison consistent, we still have 
to consider the effect of the lock-in demodulation filter 
which has a transfer function TuA^He) = 0.82, the noise 


ex 7LiA(w H e) « 35fT/Vife. We 


before the lock-in is thus pg - / lia 
find that the observed noise floor agrees with the expected 
shotnoise limit. Using Eq. Q, (101 and ( [T7| , we calcu¬ 
late the expected sensitivity parameter in the shotnoise 
limit from the above value of pg N and the 3 He-FSP ampli¬ 
tude (a = 4.15 pT), yielding p% N = a^ N {p§ N , a) ■ T^ 2 « 
111 pT • s 3 / 2 . If we compare the measured sensitivity from 
Eq. ( ]l7| ) and the shotnoise limited sensitivity we find 
i)““ ~ r/g N . We can thus state that the measurement 
was shotnoise limited. In the same way we can compare 
our results to the estimated sensitivity of [19]. We use the 
magnitude of the magnetic field produced by the 3 He FSP 
oct = 6 pT reported in their paper and p e g P = 27 fT/ a/Hz 
to find r/ ss 77 pT • s 3 / 2 . This value is comparable to the 


estimated sensitivity of H3I given at the end of Sec. |2.1 
we thus conclude that the two measurements were equally 
sensitive. We note that because the measurements de¬ 
scribed in m were done in a very weak magnetic field, 
the Larmor frequency is only u>He/27r = 3 mHz. At this 
low frequency the bandwidth limitations imposed by the 
readout magnetometer are not relevant. 

We will now estimate the ultimate sensitivity r] 7 § m , 
under the assumption of a perfectly stable magnetic field. 
This limit is reached by maximizing SNR = a rms /pg 
in Eq. (10). We thus consider a 3 He FSP with maxi¬ 
mum amplitude, measured by a shotnoise limited CsOPM 
with minimal NEM. For our best paraffin-coated cells, 
operated under optimized conditions, a minimal NEM of 
~ 7fT/v / Hz has been reported [20]. The FSP amplitude 
is maximized for 100 % 3 He polarization. We further con¬ 
sider a detection of the FSP without the bandwidth lim¬ 
itations imposed by the CsOPMs driven in the fixed fre¬ 
quency mode. This could be achieved by a CsOPM driven 
in the phase-stabilized mode by a PLL with high band¬ 
width (or by a self-oscillating Cs magnetometer [33]) ■ For 
the ~ lmbar cell used this corresponds to = 41pT at 
the CsOPM position. Combining and NEM mm from 
above leads to a maximum SNR of ss 5800 \/Hz. Using 
Eq. (10) we thus find 77 ™" ss 2 pT • s 3 / 2 for the intrinsic 
sensitivity of a combined 3 He/Cs magnetometer in this 
configuration. 


7 Conclusion 

We have described the design of a compact 3 He/Cs mag¬ 
netometer prototype for absolute measurements of weak 
magnetic fields at room temperature, and have investi¬ 
gated its performance. It was shown that the magnetome¬ 
ter is capable of performing CRLB limited measurements 
within the constraints imposed by the stability of the ap¬ 
plied magnetic field. We demonstrated that a combined 


3 He/Cs magnetometer consisting of a 3 He cell and a single 
CsOPM as readout-magnetometer can measure the abso¬ 
lute value of a 1.1 pT magnetic field with a standard un¬ 
certainty of 100 fT in a measurement time of 100 s, which 
corresponds to a relative error below 10” 7 . Measurements 
with simultaneous readout by multiple CsOPMs were pre¬ 
sented and show that the standard error of the weighted 
mean field estimate decreases statistically with the num¬ 
ber of CsOPMs, reaching AB « 60 fT in 100 seconds. This 
result is important because it implies that the magneto¬ 
metric sensitivity of a combined 3 He/Cs magnetometer 
fulfills the requirements of the current and future nEDM 
(and other) experiments by applying a suitable number of 
readout CsOPMs. A gradiometric measurement was pre¬ 
sented in which common-mode noise could be suppressed 
in the differential signal. We estimated the intrinsic sen¬ 
sitivity of the prototype magnetometer for single CsOPM 
readout, corresponding to a standard of field estimation 
error of 2fT in a measurement time of 100 s. For the 
nEDM application larger 3 He cells are foreseen which will 
also increase the FSP amplitude compared to the cell used 
here where CsOPM- and 3 He-cell were of similar size, lead¬ 
ing to an even higher sensitivity. The results show that 
combined 3 He/Cs magnetometers are suitable for the ab¬ 
solute, high precision determination of magnetic fields in 
fundamental physics experiments. 
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